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The work described in this report was performed by the DSN Operations 
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Abstract 


This report presents the development of angular tropospheric refraction models 
for optical and radio frequency usage. The models are compact analytic functions, 
finite over the entire domain of elevation angle, and accurate over large ranges 
of pressure, temperature, and relative humidity. Additionally, FORTRAN sub- 
routines for each of the models are included. 
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New Optical and Radio Frequency Angular 
Tropospheric Refraction Models for Deep 
Space Applications 


I. Introduction 

There* exists here at the Jet Propulsion Laboratory 
(JPL), and particularly within the Deep Space Network 
(DSN), a need for accurate, yet modestly sized, optical 
and radio frequency (RF) angular tropospheric refraction 
models. The basic angular refraction model (and several 
close variants) currently in use at JPL consists of three 
radically different analytic functions, each applicable 
over a different range of zenith angle (zenith angle — 90° 
— elevation angle) and is therefore immediately rather 
cumbersome. Furthermore, the accuracy of the current 
JPL refraction model is not well documented, and is thus 
subject to considerable doubt. 

The present time is particularly well suited to reexam- 
ine the question of an angular refraction model for the 
following reasons: 

(1) The remote site Antenna Pointing Subsystem (APS) 
is currently being redesigned, thus affording the 
capability to easily change the angular refraction 
modeling. 

(2) The recent advent of X-band capability, with an 
antenna beamwidth of approximately 0.020°, has 
underscored the need for high-accuracy angular 
predicts. 


The angular refraction model (or variants thereof) cur- 
rently in use at JPL is as follows: 

(1) For Z< 80.26 , 

R (NVlO)tanZ 

(2) For 90 >Z >80.26°, 



(3) For Z > 90°, 


" - (iio) (° 60974 - fr 1 } [ v - J]) 

where 

R = refraction correction, rad 
Z = zenith angle (actual), deg 
Z* = zenith angle (actual), rad 
N = “refractivity" 

To gauge tin* degree of error inherent to the current 
JPL refraction model, it has been contrasted to a con- 
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tinuous set of refraction data as computed from the work 
of B. Garfinkel (see Refs. 1 and 2), and is seen in Fig. 1. 
Hie Garfinkel data is for pressure P = 760 mm of Hg and 
temperature T = 0°C; the JPL model data has been 
matched to these conditions by setting N = 288 (i,e. t 
additionally assuming relative humidity RH = 0.0%). 
The most distressingly obvious flaws in the current JPL 
model are the discontinuities in R at the two breakpoints, 
these discontinuities (and hence errors in one or the other 
segment) amounting to approximately 30 and 300 arc 
seconds (sec), respectively. (Note: For the duration of 
this report, refraction quantities will be dealt with in 
terms of arc seconds, with 0.001° = 3.6 sec.) Further 
examination of the current model discloses that the first 
two segments are dependent upon the “refractivity" N, 
and hence pressure and temperature, while the third 
segment is not. Given that the current JPL model is in- 
accurate, has very for^r* discontinuities at the segment 
brnkpoints, and is fundamentally cumbersome because 
of the tri-segment construction, it would seem to be a 
likely candidate for a more , ccurate and reasonable 
replacement. 

The approach adopted here will be to construct first 
a very accurate optical angular refraction model of the 
form 

Hop — R op ( P ,T,Z) 

and then translate these results to the radio frequency 
case by explicitly defining a function of f such that: 

R*r = Ror(P,T,Z){f(P,T,RH)} 

Section II will generate the optical refraction model, 
while Section III will construct the RF refraction model. 


II. A New Optical Angular Refraction Model 

A. General Approach to a New Optical Angular 
Refraction Model 

In the previous section, the undesirability of the cur- 
rent JPL angular refraction model was demonstrated; in 
this section the general philosophy used to generate a 
new optical angular refraction model will be dealt with. 
One starts with the fact that angular refraction is crucially 
important in the effectuation of various astronomical en- 
deavors, and hence there exists copious amounts of refrac- 
tion data. The main drawback to these data is, however, 
that they are either in tabular form or are calculated via 
schemes which require large amounts of tabular inputs 


(for instance, see Refs. 1, 2, 3, and i). Furthermore, the 
astronomical accuracy requirements are very stringent 
(perhaps down to about the 1-sec level), while the DSN 
requirements ire no greater than the 10-sec level. In light 
of the above, it is clear that one reasonable approach 
would be to use empirical methods to develop a simple 
analytical expression to approximate the very accurate 
astronomical refraction data available. Since the envi- 
sioned use of the new model within the DSN includes 
small i emote site computers as well as large central com- 
plex computers, desirable features would include: 

(1) A single expression over the entire domain of Z, 
instead of multiple segments, each applicable over 
different ranges of Z. 

(2) Accuracy to about the 10-sec level for reasonable 
ranges of Z and tropospheric conditions. 

(3) Model to be designed to minimize both computer 
memory and run time. 

B. Selection of an Optical Angular Refraction 
Data Base 

After a review of the literature, it became apparent that 
a reasonable selection for a data base would be the work 
of Boris Garfinkel of the Yale University Observatory 
(see Refs. 1 and 2). Garfinkel’s original theory was pub- 
lished in 1944, and then reexamined in 1966. The form of 
his model is semi-analytical in that it is a closed function 
with Z, P, T as variables, but also requires tabular input 
in the form of Z-dependent polynomial coefficients. 
More importantly, his model is continued for Z > 90°, an 
aspect which is most frequently missing. in Other angular 
refraction works. Finally his work compares well with 
other authorities in the field. For instance, Garfinkel com- 
pares his data at P — 760 mm and T = 0°C with those of 
the Radau and the Pulkova models as follows (with 2 V 
(observed zenith angle) in degrees and R in seconds): 


Z' 

R 

Garfinkel 

R 

Radau 

R 

Pulkova 

80 

331 

331 

331 

81 

,368 

366 

385 

82 

410 

409 

408 

83 

463 

462 

460 

84 

531 

529 

527 

85 

619 

617 

614 

86 

738 

735 

733 

87 

905 

903 

900 

88 

1153 

1152 

1147 

89 

1544 

1545 

1537 

90 

2206 

2208 

2199 
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Table 1 provides a detailed tabulation of Garfinkel 
refraction data for 0° < Z < 93°, P = 760 mm, and 
T = 0°C. 


where 


K| — 7T [(Ri)ma* "b (fli)mln] 


C. Derivation of a Basic Optical Angular 
Refraction Model 



(fl.)u 


The needs of the DSN for an angular refraction model 
are restricted to the following range of Z': 



K.l 


where 


0° < Z' < 92° 


and where the conditions satisfied are the following 
n + 1 equations in n + 1 unknowns. 


Z' — Z — R(Z) = observed zenith angle 

this range being encompassed by the Garfinkel data in 
Table 1. It was hoped that the data base chosen (i.e., the 
Garfinkel data whose selection was discussed in the pre- 
vious Subsection) could be approximately fit to a function 
(or functions, as necessary) via routine least squares tech- 
nique's. It was planned to do all work for P = 760 mm 
and T = 0°C under the assumption that P and T effects 
could Ik* (multiplicatively) added at a subsequent time. 
The data base chosen was a slightly smaller subset of the 
data base displayed in Table 1. The frequency of data 
points was rather arbitrarily chosen as follows: 


Let 


A, = R,( Z.) - X(Z.) 

.-[&«]* 

Then, finally. 


da 


= o 


Pa 

8* 


= 0 


Range of Z, deg Data frequency, deg 


0 < Z < 70 

0.5 

70 < Z < &5 

0.2 

85 < Z < 93 

0.1 


with the net effect that the refraction data were increas- 
ingly "weighted” in the high Z region where the rate of 
change of refraction is the greatest. The computer pro- 
gram utilized in this study is a standardized least squares 
subroutine available to all UN1VAC 1108 users at JPL 
(see Ref. 5). Basically, it fits a data set to an nth degree 
polynomial such that the residuals are minimized in a 
least squares sense, i.e., 

if R,(Z,); i - data set 
then a function X is fnrmulah'd such that 


It was originally intended to attempt a least squares 
curve fit to the "raw” refraction data, shown in Fig. 2. It 
was observed, however, that the natural log (In) of R gave 
a very smooth representation and possessed, of course, 
far less dynamic range, as can Ik- seen in Fig. 3. It seemed 
possible that it might yield a Ix-ttcr fit for a lower order 
polynomial (a desirable property), i.e.. fitting: 

In (R,(Z,)); X 

Finallv, it was observed that taking the Inverse tangent 
(arctan) of In (R) yielded a representation that appeared 
almost linear, as can be seen in Fig. 4. This was also felt 
to be worth attempting as a fit, in the form of: 


X = ±K h ,{V( Z))> 

fmO 


I In (fl,(Z,)) 

I In (R(45 0 ))/ 


X 
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The* In (R) fit was attempted first as the most likely 
candidate. The main goal established was to find the 
smalle'st order fit which would ke*e*p the maximum resid- 
ual below some reasonable limit. The results of a se*ries 
of different orde r fits appear in Fig. 5. Although it might 
at first serin strange that the* absolute maximum residual 
does not decrease monotonicallv with ele*grer of fit. all one 
should really expe-ct is that <r decrease monotonicallv with 
de*gre*e* of fit, and this was the* case. At anv rate, the* 8th 
degree case* was fe*lt to lx* the* best ceimpromise, as one 
needed to gei to a 14th ele-grer fit to obtain significant im- 
provement. The* e>the*r two types of fits were* atte*mpte*d, 
with a general tradeoff e*\perte*el that an increase in func- 
tional depe*nde*nce (In, arctan, etc.) should decrease the 
order of fit necessary. A few of the* salient fe*ature*s of each 
of the* thre*e* tvpes of fits attempted are: 

(1) Raw Data Fit 

(a) Simple* — no additional functions required for 
modeling. 

(b) Minimum aeceptable polynomial required 
~12th order 

(c) l,urge residuals ( — 20 sec eir higher) as Z— *0°, 
leading te> unpalatable* result of re*fraction lx*ing 
applied in wrong direction at very small Z, e*tc. 

(2) Ln(fl) Fit 

(a) Model weiulel re*e|uire* natural oxjxmentiation 
(exp). 

(b) Minimum acceptable polynomial required 
—8th order. 

(c) logarithmic condition of fit forces residuals to 
be* approximately proportional to R. so resid- 
uals ejuite* small exce*pt at very large* Z. 

(3) Arctan (ln(/?)) Fit 

(a) Model would require tangent (tan) and exp. 

(b) Minimum polynomial fit re*ejuire*d —6th order. 

(c) Extremely low residuals for Z '-c 90° and ejuite 
high residuals for Z ^ 90 . 

A conqiarison of the tlire*e* types of fits is seen in Fig. 6.' 
The In R fit was assessed to lx* the* best compromise* 


'This figure and Kilt*. 7. 9. and 10 wore prepared on the basis of 
interim results and are at variance with the final model liv as 
much as 5 sec at large Z. Therefore, they should lie used for 
illustration only. 


amongst the design goals stated in Subsection A. Further 
refinement to the 8th order In R fit was accomplished by 
making minor adjustments to the data set used in the fit 
process, until an optimum fit (in the sense of the smallest 
maximum residual) was achieved. For mis case, the maxi- 
mum residual in the interval 0° < Z < 92° occurred at 
about Z = 91.1° and had a value of: 

A R - +21.6 sec 

D. Complete Optical Refraction Model Determination 

The refraction model, as finally determined in the pre- 
vious section, is as follows: 

n = ex P {tK,.,[I7(Z)]4-K„ 

where 

R = refraction, sec 
Z = zenith angle, actual 
EL = 90° — Z = elevation angle 

K x = 46.625 
K, = 45.375 
K, = 4.1572 
K, = 1.4468 
K, = 0.25391 
K, = 2.2716 
K t = -1.3465 
K, = -4.3877 
K v = 3.1484 
Km = 4.5201 
K,i = -1.8982 
= 0.89000 

When this model is compared to the (•arfinkcl data (with 
/’ 760 nun and T 0 C), the following maximum 

residuals* rr*sult: 

0 <Z<85° Afl = + 5.6 sec 

•'All residuals i A/t i will b<* (>irfinkd D.it.i — l*ro|Niwd Model 


4 
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85° <Z<92° aR = + 21.8 sec 

92° < Z < 93° aR = —302.6 sec 

The very large residuals between Z = 92° and Z = 93° 
are primarily a result of ending the fit at 92°. At this 
point there was still one point of concern and that was: 

As Z > 93° 

| R | — * very largo 

such that at some Z>93® there exists the following 
condition: 

Z - R(Z) < 90° (or local horizon) 

giving the appearanc*- of a “false” rise. This would, of 
course, jxise difficulties for trajectory-type programs which 
calculate, and examine for a rise condition, zenith angles 
considerably larger than 90°. Based on this undesirable 
feature, it was felt that the model should be modified 
such that shortly after Z = 93° it would Ik* required that: 

RfZ)-»0 

At the same time, it was felt that possibly the charac- 
teristics of the model for Z ^ 90° could be improved 
upon. The residuals between Z 90° and Z = 93° look 



It was therefore felt that if a function (sav, A,(Z)) could 
lie derived with inverse characteristics to tin above resid- 
uals plus |M>ssessing the following qualities: 

A (Z) — * very’ large for Z j 93 

A ,(Z) — » very small for Z V. 9) 

then a model of the form: 

fl - exp [IW) /[l + A,(Z)]J K ia 


could perhaps both improve the present model between 
90 c and 93° and drive the model to approximately zero 
(actually 1 - K,,) thereafter. A function to accomplish 
this was constructed (empirically) as follows: 

Aj(Z) = (Z — C«) {exp [C,(Z — C»)]} 

where 

Z = zenith angle, deg 
Co = 91.870 
C, = 0.80000 
C, = 99,344 

Tlie improvement in the Z = 90° to Z = 93° region 
can be seen in Fig. 7, while the rapid drop off of the 
modified model after Z 93° can be viewed in Fig. 8. 
The maximum ri-siduals after the above modification 
become: 

0° < Z < 85“ aR = + 5.8 sec 

85° < Z < 92° AR = - 14.7 sec 

92 <Z<93° aR ~ 15.0 sec 

E. Refraction Model Functional Dependence 

Upon Pressure, Temperature, and 

Relative Hun ddity 

It was originally felt that once a refraction model for 
standard conditions (P 760 mm and T - 0 e C) had been 
achieved, the usual scaling bv /’/760 and 273 A T 4 273) 
could be applit-d However, after examining differenr 
combinations of P and T in the Garfinkel data, this did 
not prove to be an adequate treatment of the pressure 
and tein|M-rature dependence, and additional work in this 
area was required. 

1. Pressure correction. Examination of the Garfinkel 
data at different pressures indicated that scaling of the 
basic model by PI 700 was reasonable at most Z. but 
broke down as Z ^ 90°. It was Imped that this could be 
compensated for by a correction factor (say, A ,(P. Z)) such 
that the entire pressure correction factor would 1 m- of the 
form: 

4j<‘ -•M/’.Z)} 
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Furthermore, it would be necessary that: 


in combination with the conditions: 


A,(P,Z)-*0. Z<90° 

A,(P,Z)-0, Z>93° 

It was noted in the examination of the Garfinkel data 
that the pressure effect (as different from P/760) was for 
the most part separable, i.e.: 

A,(P.Z)~A,(P)A,(Z) 
and it could be seen that further: 

A,(P) ~ (P - 760) 

A representation for A.- was then empirically constructed 
as follows: 

A,(Z) ~ exp [A t (Z - A,)] 
so that the- A,i P. Z) pressure correction would be: 

Ai(P , Z) = (P - 760) exp [ A,(Z - A,)] 

The results of using A,(P, Z), above, can be seen in 
Fig. 9. Finally, to satisfy the conditions of a small A,(P, Z) 
for Z ^93 , the previously determined A,(Z) was utilized 
to arrive at the following expression: 

Pressure correction factor = 

_Pj _ (P — 760) exp [A,(Z — A a )] ) 
760} 1+A,(Z) / 

where 

Z = zenith angle, deg 
P - pressure. mm of Hg 
A, = 0.40816 
A, = 112.30 

A|(Z) = as previously defined 

2. Temperature correction. The investigation of tem- 
perature effects proceeded along the same lines as the 
investigation of pressure effects in the previous section, 
with the goal of a total temperature correction factor 
in the form of: 

273 

T ‘ 273 {1 “ 


x(T.Z)-*o z< 90 ° 

A,(7\Z)-»0, Z>93° 

Similarly, the temperature effect was found to be ap- 
proximately separable: 

A,(7\Z)~A,(T)A.-(Z) 

and the following was (empirically, constructed: 

A ,~T 

X ~exp[B,(Z - B,)] 

so that the A ,(T. Z) temperature correction would be: 

A,(T, Z) = (T) exp [B,(Z — B,)] 

The results of using A,(7\ Z), above, are seen in Fig. 10. 
Once again, to satisfy the condition of a small A,(T, Z) 
for Z > 93°, the previously determined A,(Z) is utilized 
to arrive at the following total expansion : 

Temperature correction factor = 

273 f (T)exp[B,(Z- B,)» 

7^ 2731 1 + A,(Z) / 

where 

Z = zenith angle, deg 
T = temperature, ®C 
B, = 0.12820 
B, = 142.88 

A,(Z) = as previously determined 

3. Relative humidity correction. Both Garfinkel (see 
Ref 2) and the Pulkova Model (see Kef. 4 i.idicate that 
the correction for relative humidity is veiy small, per- 
haps on the order of several seconds at large Z. at a 
maximum It was therefore decided here to ignore cor- 
rections based on relative humidity, at least until some 
time when a stronger ease can Ik- math- for the necessity 
of its inclusion, given the level of accuracy (10 six) 
inherent to the model here ticing proposed. 
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F. Complete Optical Angular Retraction With Pressure 
and Temperature Corrections 

The final refraction modrl with pressure and tempera- 
ture* accounted for is as follows: 


R = F p f\cx p)> ! + A l( Z) ) 




z) i\ 
i + a,(z )/; 

X(T,Z) |\ 
1 ' 


A,(P.Z) = (P - P 0 ){exp [A,(Z — A,)]} 
a,( 7\Z) = (T — '** 0 ){exp [B,(Z - Bj)]) 


A, = 0.40810 

A, = 112.30 

T = temperature, K 

T„ ■■= 273.00 K 

B, = 0.12820 

B, = 142.88 
Co = 91.870 

C, = 0.80000 
r. = 99.344 


The accuracy of this model for various pressures, tem- 
perature, and ranges of Z. as compared to the Garfinkel 
data, can be seen in Table 2. 


Aj(Z) =(z-c„v: *xp[c,(z c 3 )]i 


where 

R — refraction, sec 
Z actual zenith angle, deg 
EL = 90 Z elevation angle 


The signature of the residuals at large 7 and with 
P 760 mm and T 0°C can be seen in Fig. 7 (modified 
8th order In R fit). For use where simplicity is of a more 
urgent need than accuracy, an abbreviated version of the 
model can be obtained l>\ setting: 

A, = A- = A, = 0 



*■ < 
Kj 1 


AC, = 46625 
K, = 45.375 
K, = 4.1572 
AC, 1 4468 
AC- = 0.2539! 
K. = 2 2716 


such that one has: 

where all quantities are as previous!) defined. The accu- 
racy of this abbreviated version, once again as compared 
to the Garfinkel data, is seen in Table 3. Also, the effects 
jf the deletion of A is seen in Fig 9, of A, in Fig. 10, and 
of A , in Fig. 7. 


K 1.3465 

K. = 4 3877 

K, = 3.1484 
AC,.. = 4.5201 
AC,, 1.8982 


Q. FORTRAN Subroutines of the Optical Refraction 
Models 

Appendix A pri-si nts a FORTRAN subroutine of the 
full mode) described in Subsection F, while Appendix B 
presents a FORTRAN subroutine corresponding t<> tin 
abbreviated model, also describ'd in Subsection F In 
puts required are as follows: 


K,, - 0.89000 PRESS = pressure, nun of Mg 

P — pressure, mm of Mg TEMP = temperature, K 

P„ = 700.00 mm ZN1TH — actual zenith ungle. deg 
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and the subroutines) return with: 

R = refraction correction, sec 

III. A New Radio Frequency Angular Refraction 
Model 

In the previous section, a new optical angular refrac- 
tion model was presented. In this section, past attempts 
to translate from an optical model to a radio frequency 
model will he dealt with, and then a new method to 
accomplish this transformation will be proposed. 

A. Past Attempts to Transform Angular 
Refraction Models From Optical to 
Radio Frequencies 

To facilitate a discussion of past attempts to generate 
angular tropospheric refraction models for use at radio 
frequencies, let the following notation be introduced: 

Rup = R(P,T,Z) = optical refraction model from 
section II 

R hr = Rnr(l\T.Z,RH) = radio frequency refrac- 
tion model 


this expression by the total radio frequency surface 
refractivity, i.e., 

refraction ~ 

where N* = reference optical refractivity. 

At this point, one must ask, what are the implications 
of this procedure? Since any signal (that is of interest 
here) must traverse the entire troposphere, and is of 
course, continually being refracted, one might think that 
instead of being proportional to surface refractivity, angu- 
lar refraction is really more nearly proportional to total 
(integrated) tropospheric refractivity, i.e., 



However, refraction could also be proportional to surface 
refractivity if it could be assunml that there exists some 
f(h) such that: 

N(/.) =; N.f(h) 

so that 


P — pressure 
T = temperature 
Z = zenith angle 
RH — relative humidity 
N(h) = ND(h) 4 N\V(h) 

N(h ) = total refractivity at radio frequencies 
ND(h) = dry, or optical component, of refractivity 
NW(h) = wet component of refractivity 
h = height 
h 0 = station height 


refraction oc f N,f(h) dh — N, / f(h) dh 

Making the assumption that 

ND(h) ~ ND./,(h) 

NW(h) ~ NW.Hh) 

one would have for the optical case: 

R, oc / ND(h) dh = J ND.fdh) dh 
= ND.ff x (h)dh. 

For the radio frequency case: 


s = parameter surface value 

ND(h 0 ) = ND. 

NW(h n ) = XW. 

N(h 0 ) = N. 

In general, attempts to construct a radio frequency 
refraction model consisted of appropriating an empirical 
model from optical refraction work which would give the 
functional dependence on Z (say R,( Z)), and then scaling 


/W * fN(h)dh - / (XD(h) 4 NW(h)) dh 

= / {ND.f,(h) 4 NW.f,(h)) dh 
= ND, f f,(h) dh + N\V, / f .(h) dh 

Without precise knowledge of tin* form of f x (h) and /.(/»), 
the only way that the surface refractivities could be used 
to transform from tin- optical case to the radio case would 
be if 

f,(h)~h(h) 
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Then 

R„,« (ND M + NW,)ffi(h)dh 

and indeed 

R„, oc(ND, + NW.)[-^] 

However, it is well known that wet refractivity “decays” 
much more rapidly than dry refractivity (for instance, 
Ref. 6), so that f,(h) and f 2 (h) are quite dissimilar; thus, 
the procedure of scaling an optical angular refraction 
model by the total surface radio refractivity to achieve a 
radio angular refraction model would appear to be seri- 
ously flawed. 

B. Method Used to Transform From an Optical 
to a Radio Frequency Refraction Model 

From the previous section it was seen that 
Rkp « N t = ND„ + N\V, 
is a poor choice. A more logical choice would be 

JW « SN(h) dh = J [ND(h) + NW{h)} dh 

— f ND(h) dh + / NW(h) dh 

f. . SNW(h)dh ) 

- / nd ( h) *|1+ 

Similarly, for the optical case (using the model previously 
presented) : 

R 0P oc J ND(h) (Ui 

Combining the above, one arrives at the equation thet 
will be used for the radio frequency angular refraction 
model: 

ft„{P,7\Z,RH) = R„,(P,T,Z) j 1 + j^hldh } 

C. Determination of Ratio of Integrated 
Wet Refractivity to Integrated Dry 
Refractivity 

In attempting to determine an analytical parametric 
representation for the expression: 

fNW(h)dh 
/ ND(h) dh 


the most difficult problem by far lies with the integrated 
wet refractivity. Berman first showed in 1970 (Ref.7) that 

/ND(h)dh = AP,[|-] 

where 

A = 77.6 

P, = surface pressure, mbar 
R = perfect gas constant 
g = gravitational acceleration 
g/R = 34.1°C/km 

and also gave an expression to approximate the inte- 
grated wet refractivity: 

fNW(h)dh = 



where 

Cx = 77.6 
C 2 = 29341.0 
RH = relative humidity 
y ~ temperature lapse rate 
C = 38.45 

T 0 — extrapolated surface temperature 
A = 7.4475 In (10) 

B = 2034.28 In (10) 

Chao (Ref. 8) later improved upon the integrated wet 
refractivity with the expression: 

f NW(h) dh = 1.63 X 10- + 2.05 X 10 2 « 
where 

e 0 = surface vapor pressure, N/m* 

T„ = surface temperature, K 
a = temperature lapse rate, K/km 
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However, both of these expressions depend upon one or 
more parameters not measurable at the surface (i.e., 
temperature lapse rate, etc.), and neither is particularly 
accurate, doing back to the previous section, if the 
altitude-dependent refractivities could really be repre- 
sented as 


XD(h) - ND.f t (h) 



h < 12.2 km 

iP)} 

h < 13 km 


NW(h) - NW.f a (h) 


— 0 h > 13 km 


and if the above refractivities eow.d Ik* integrated, i.e., if 
A below could be evaluated as 


Performing the dry rcfractivity integration, one has 


, _ SMh)dh 

JfJFuIK 

then one might simply expect that 

/ N\V(h) dh iNW.) 
f ND(h) dh ~ A \ ND,f 

To test this hypothesis, the authors had ten cases used in 
Ref. 7. Although a very small number, the cases were 
alternate day and night profill’s selected throughout the 
year (December, February, April, August, September). A 
least-squares linear curve fit to the above data was per- 
formed as follows: 


+ /* m m - i exp (" el’ 21 )} ih 


+ 


70 xrrs Z -50 r (h — 12.2) 1 

-gj— \< lh 

transforming the first integral by 


fNW(h)dh jNW.) 

jXD(h)dh ’ A \ ND, f 

Tlie fit yielded the following: 


A = 0.3224 
,t{%) = 00.93% 

/ ( NW(h)dh A ( NW. 1 \ 

<t(%) - ioo > * y (ND(h)( i h - A | VD ( j 


Translated to centimeters of integrated rcfractivity, one 
would have 


fr(cm) — 2.0 cm 

Table 4 and Fig. 1 1 present the detailed analysis of the 
ten cases described. 




dh = — 42.7 dx 


so that 


/( 1 - 4 J 7 ) Hh = -42.7/Vifc 


= - 42 - 7 tt 




= 6.952 


Transforming the second integral by 


As a totally independent check of this observed rela- _ ]2 2) 

tionship, use can be made of work done by Chao (Ref. 6) ^ = x 

on wet and dry rcfractivity profiles. Combining Eqs. (9). 

(10), (13), (14), (15). and (16) from Ref. 6 one has: dh = — 6.4 dx 
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so that 

Combining the integrated wet refractivity and the inte 
grated dry refractivity yields 

/ exp [~ 6 4 ] dh = 64 / ex P ( x ) dx 

= — 6.4 exp (x) 

f„ NW W d*' 2.6 (NW.) 

‘ a6m(ND - ) 


» CS) 

= 6.4 



This is to he compared to the previously determined 
Or, finally relationship from actual data of: 


/: 


ND(h)dh=ND,( 0.952) 


70 

+ 209 * 0 .( 6 . 4 ) 


= 8.6174(*VD„) 


Performing the wet refractivity integration, one has 


/; 


NW(h) dh = 


/'W.(i-A)',". 


f NW(h)dh 

Jji ~ n <ic 

r 


= 0.3224 


ND(h) dh 



Since the value of the l.r standard deviation 


l«r - 00.93% (-2 cm) 

found from actual data compares favorably with the most 
recent modeling published by Chao in Ref. 8 (~3 cm for 
combined night and day profiles), and since the basic rela- 
tionship seems verifiable by average profiles presented by 
Chao; the determined expression will be adopted for use 
with the optical refraction model. Tin- surface refractivity 
( Ref. 7) is defined as: 


Transforming the integral by 

(>-*)-■ 

dh = -13 dx 

/( 1 -ySj) dh = - 13 Jx 1 dx 


so that 


= - 13 t 




= 2.6 


J* NW(h) dh = 2.6 (NW.) 


NW. = 


(RH).CtC. 


T: 


exp 


(AT. B\ 
\T. Cj 


ND. = C, ~ 
so tint one would obtain 


j, +(o ^^ 

| I ND(h) dh | * * ' ’ 

To integrate this expression into the optical model, the 
pressure term must be converted from mbar to mm: 

P.(mbar) = P.(rnm) X -=^r 

(Oil 

so that one would finally have 

|, / NW(h) (lh\ , , (7.1X10") (RH). (AT. B\ 

\ l + jrnmar 1 + — tt. — , ' xi> (fr - r) 
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where 

( RH ), = surface relative humidity (100% = A .0) 

T, = surface temperature, K 
P, = surface pressure, mm of Hg 
A = 17.149 
B = 4684.1 
C = 38.450 

D. Final Angular Tropospheric Radio 
Frequency Refraction Model 

The following gives the complete radio frequency 
angular tropospheric refraction model: 


W £k,. 5 [[/(Z)]>) \ 

' "i + a j( z) i - *■■/ 


F P = 
F,= 
F„ = 


(p_( a,(p,z) n 

\p.. I 1 + A,(Z )f) 

(hU- ^ X) 

\T\ 1 * A ,(Z) f ) 

(. W 0 RH ( rwj 

v + ~TP~~\ exp T 3 


-W, 


W 3 



A|(P,Z) = (P - P n ) {exp [A,(Z — A : )]} 
A 2 (T,Z) = (T-T„){exp[B 1 (Z-B : .)]} 
A,(Z) = (Z - C„) {exp [C,(Z - C a )]> 
where 

R = refraction, sec 
Z = actual zenith angle, deg 
EL = elevation angle 
EL = 90 deg - Z 



K, = 46.625 
K 2 = 45.375 
K, = 4.1572 
K, = 1.4468 
K> = 0.25391 


Ko = 2.2716 
K ? = — 1 .3465 
K t = -4.3877 
K, = 3.1484 
K 10 = 4.5201 
= -1.8982 
= 0.89000 

P = pressure, mm Hg 
P 0 = 760.00 mm Hg 

A, = 0.40816 
A* = 112,30 

T = temperature, K 
T u = 273.00 K 

B, = 0.12820 
B 2 = 142.88 
C„ = 91.870 

C, = 0.80000 
C 2 = 99,344 

RH = Relative humidity (100% = 1.0) 

W 0 = 7.1 X 10' 

W, = 17.149 
W, = 4684.1 
W, = 38.450 

E. Model Accuracies 

Tin- inaccuracies introduced by the wet rcfractivity 
term predominate over the inaccuracies presented in 
Section II. Considering 

l<r = 1.00% 

the maximum l<r angular errors would be 


Z. deg 

Afl, sec 

AR, deg 

0-85 

6 

0.002 

85-90 

18 

0005 

90-93 

50 

0.015 
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F. FORTRAN Subroutines of the Radio Frequency 
Refraction Models 

Section II presented two FORTRAN subroutines, cor- 
responding to the full optical refraction model and an 
abbreviated version. These two routines have l>een trans- 
formed to the radio frequency version of the refraction 
model, and are presented in Appendixes C and D. The 
FORTRAN subroutine SBEND (Appendix C) represents 
the full model, while XBEND (Appendix D) gives the 
abbreviated version. Inputs required are: 


PRESS = pressure, mm of Hg 
TEMP = temperature, K 
HUMID = % of relative humidity (100% = 1.0) 
ZNITH = actual zenith angle, deg 

and the subroutines return with 

R = refraction correction, sec 
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Table 1. Garfinkal refraction data* for P = 760 mm and T = 0*C 


z 

R 

R' 

Z 

R 

R' 

Z 

R 

R' 

0.0 

0.00 

0.00 

23.5 

26.33 

26.34 

47.0 

64.95 

64.99 

0.5 

0.55 

0.55 

24.0 

28.97 

28.97 

47.5 

66.10 

68.15 

1.0 

1.07 

1.08 

24.5 

27.60 

27.81 

48.0 

67.28 

87.32 

1.5 

1.61 

1.61 

25.0 

28.25 

28.26 

48.5 

68.47 

68.52 

2.0 

2.15 

2.15 

25.5 

28.90 

28.91 

49.0 

69.89 

69.73 

2.5 

2.68 

2.68 

26.0 

29.55 

29.58 

49.5 

70.93 

70.97 

3.0 

3.22 

3.22 

26.5 

30.21 

30.22 

50.0 

72.19 

72.24 

3.5 

3.75 

3.76 

27.0 

30.88 

30.89 

50.5 

73.48 

73.53 

4.0 

4.29 

4.29 

27.5 

31.55 

31.56 

51.0 

74.79 

74.85 

4.5 

4.83 

4.83 

28.0 

32.23 

32.24 

51.5 

76.14 

78.19 

5.0 

5.36 

5.36 

28.5 

32.91 

32.92 

52.0 

77.51 

77.57 

5.5 

5.90 

5.90 

29.0 

33.60 

33.62 

52.5 

78.91 

78.97 

6.0 

6.43 

6.44 

29.5 

34.30 

34.32 

53.0 

80.34 

80.40 

6.5 

6.97 

6.97 

30.0 

35.01 

35.02 

53.5 

81.80 

81.87 

7.0 

7.51 

7.51 

30.5 

35.72 

35.73 

54.0 

83.30 

83.37 

7.5 

8.05 

8.05 

31.0 

36.43 

36.45 

54.5 

84.83 

84.91 

8.0 

8.59 

8.59 

31.5 

37.16 

37.17 

55.0 

86.40 

86.48 

8.5 

9.12 

9.13 

32.0 

37.89 

37.91 

55.5 

88.02 

88.10 

9.0 

9.67 

9.67 

32.5 

38.63 

38.65 

56.0 

89.68 

89.77 

9.5 

10.21 

10.21 

33.0 

39.38 

39.40 

. 56.5 

91.39 

91.48 

10.0 

10.75 

10.75 

33.5 

40.14 

40.16 

57.0 

93.14 

93.23 

10.5 

11.30 

11.30 

34.0 

40.90 

40.92 

57.5 

94.94 

95.03 

11.0 

11.84 

11.84 

34.5 

41.68 

41.70 

58.0 

96.78 

96.88 

11.5 

12.39 

12.39 

35.0 

42.46 

42.48 

58.5 

98.67 

98.78 

12.0 

12.94 

12.94 

35.5 

43.26 

43.27 

59.0 

100.61 

100.72 

12.5 

13.49 

13.49 

36.0 

44.06 

44.08 

59.5 

102.01 

102.72 

13.0 

14.04 

14.05 

36.5 

44.87 

44.89 

60.0 

104.66 

104.78 

13.5 

14.60 

14.60 

37.0 

45.69 

45.71 

60.5 

106.76 

106.89 

14.0 

15.15 

15.16 

37.5 

46.53 

46.55 

61.0 

108.93 

109.00 

14.5 

15.71 

15.72 

38.0 

47.37 

47.39 

81.5 

111.17 

111.31 

15.0 

16.28 

16.28 

38.5 

48.23 

48.25 

62.0 

113.49 

113.64 

15.5 

16.84 

18.85 

39.0 

49.10 

49.12 

82.5 

115.93 

116.10 

16.0 

17.41 

17.41 

39.5 

49.98 

50.00 

63.0 

118.47 

1 18.04 

16.5 

17.98 

17.98 

40.0 

. 50.87 

50.89 

63.5 

121.10 

121.28 

17.0 

18.55 

18.56 

40.5 

51.77 

51.80 

64.0 

123.79 

123.98 

17.5 

19.13 

19.13 

41.0 

52.09 

52.72 

64.5 

126.53 

128.72 

18.0 

19.70 

19.71 

41.5 

53.62 

53.65 

65.0 

129.35 

129.56 

18.5 

20.29 

20.29 

42.0 

54.56 

54.59 

65.5 

132.28 

132.50 

19.0 

20.87 

20.88 

42.5 

55.52 

55.55 

68.0 

135.33 

135.57 

19.5 

21.46 

21.47 

43.0 

56.50 

56.53 

66.5 

138.51 

1 . 38.70 

20.0 

22.06 

22.06 

43.5 

57.49 

57.52 

67.0 

141.82 

142.09 

20.5 

22.05 

22.66 

44.0 

58.50 

58.53 

07.5 

145.20 

145.55 

21.0 

23.26 

23.66 

44.5 

59.52 

59.56 

08.0 

148.86 

149.17 

21.5 

23.86 

23.87 

45.0 

00.58 

60.60 

68.5 

152.04 

152.97 

22.0 

24.47 

24.48 

45.5 

81.83 

61.07 

09.0 

156.01 

156.97 

22.5 

25.09 

25.09 

46.0 

62.72 

02.76 

09.5 

160.75 

161.13 

23.0 

25.71 

25.72 

46.5 

63.83 

63.87 

70.0 

105.00 

105.40 

•H gives tlu- 

retraction correction if Z 

actual while 

R' gives the refraction correction if Z 

observed. 
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Table 1 (contd) 


z 

K 

fi ' 

Z 

R 

fl ’ 

Z 

R 

R' 

70.2 

lbfl .83 

167.25 

79.8 

323.92 

328.68 

87.2 

895.71 

947.24 

70.4 

168.83 

169.06 

80.0 

330.09 

333.02 

87.3 

914.73 

969.38 

70.6 

170.47 

170.91 

80.2 

336.46 

339,52 

87.4 

934.45 

992.51 

70.8 

172.33 

172.78 

80.4 

342.99 

348.20 

87.5 

954.84 

1016.58 

71.0 

174.23 

174.69 

80.6 

349.75 

353.18 

87.8 

975.96 

1041.62 

71.2 

178.18 

176.64 

80.8 

3 . 56.81 

380.45 

87.7 

997.91 

1067.68 

71.4 

178.13 

178.63 

81.0 

364.17 

368.05 

87.8 

1020.64 

1095.01 

71.6 

180.15 

180.66 

81.2 

371.82 

375.93 

87.9 

1044. ’.6 

1123,53 

71.8 

182.20 

182.73 

81.4 

379.76 

384.12 

88.0 

1068 . S 3 

1153 28 

72.0 

184.30 

184.85 

81.6 

387.99 

392.61 

88.1 

1093.93 

1 184.47 

72.2 

186.44 

187.01 

81.8 

398.55 

401.48 

88.2 

1120.27 

1217.10 

72.4 

188.64 

189.23 

82.0 

405.48 

410.74 

88,3 

1147.59 

1251.29 

72.8 

190.90 

191.51 

82.2 

414.76 

420,35 

88.4 

1176.01 

1287.15 

72.8 

193.21 

193.84 

82.4 

424.44 

430.45 

88.5 

1205,55 

1324.68 

73.0 

195.57 

196.22 

82.6 

434,59 

441.05 

88.8 

1236.25 

1364.16 

73.2 

197.97 

198.64 

82.8 

445.21 

452.12 

88.7 

1268.19 

1405,55 

73.4 

200.41 

201.10 

83.0 

456,30 

483.71 

88.8 

1301.38 

1449.01 

73.6 

202.90 

203.82 

83.2 

467.87 

475.82 

88.9 

1335.90 

1494.78 

73.8 

205.45 

206.19 

63.4 

480 . 0.3 

488.71 

89.0 

1371.84 

1543.13 

74.0 

208.07 

208.84 

83.6 

492.90 

502.21 

89.1 

1409.18 

1594.01 

74.2 

210.75 

211.56 

83.8 

506,32 

516.28 

89.2 

1448.01 

1647.64 

74.4 

213.51 

214.35 

84.0 

520,31 

531.10 

89.3 

1488.47 

1704.25 

74.6 

216.34 

217.20 

, 84.2 

535.04 

546.76 

89.4 

1530.70 

1704.12 

74.8 

219.23 

220.12 

84.4 

550.57 

* 563.28 

89,5 

1574.641 

1827.44 

75.0 

222.18 

223.11 

84.6 

566.91 

580.78 

89.6 

1620.40 

1894.34 

75.2 

225.21 

226.18 

84.8 

584.18 

599.35 

89.7 

1668.02 

1965.25 

75.4 

228.33 

229.34 

85.0 

602.50 

619.11 

89.8 

1717.65 

2040,56 

75.6 

231.54 

232.80 

85.1 

612.07 

629.38 

89.9 

1769,36 

2123.12 

75.8 

234.84 

235.94 

85.2 

621.88 

639.99 

90.0 

1823.24 

2205,54 

76.0 

238.20 

239.34 

85.3 

631.94 

650.89 

90.1 

1879.28 

2298.34 

78.2 

241.64 

242.81 

85.4 

642,32 

662.10 

90.2 

1937.63 

2392.18 

78.4 

245.15 

246,37 

85,5 

652.95 

673.70 

90.3 

1998.35 

2495.00 

76.8 

248.77 

250.05 

85.6 

663.88 

685.73 

90.4 

2062.49 

2604.75 

76.8 

252.50 

253.85 

85.7 

675.18 

698.15 

90 5 

2130.07 

2722.08 

77.0 

256.35 

257.75 

85.8 

686.86 

7 10.99 

90.6 

2196.81 

2847,58 

77.2 

200.29 

281.74 

85.9 

898.89 

724,30 

90.7 

2269.53 

2982.00 

77.4 

284.33 

265.85 

86.0 

711.31 

738.00 

90.8 

2343.68 

3126.84 

77.6 

268.50 

270.10 

86.1 

724.15 

752.11 

90.9 

2419.93 

3282.09 

77.8 

272.80 

274.48 

86.2 

737.33 

766.87 

91.0 

2500.7 1 

3450.46 

78.0 

277.24 

279.00 

86,3 

7 . 50.87 

782.12 

91.1 

2584,52 

. 3632.11 

78.2 

281.80 

283.63 

86.4 

764.99 

797.78 

91.2 

2671.66 

3827.51 

78.4 

288.49 

288.41 

86.5 

779,57 

814.09 

91.3 

2762.19 

4039.04 

78.6 

291.32 

293,34 

86.6 

794,50 

631.04 

91.4 

2856.17 

4209.13 

78.8 

296.33 

298.45 

86.7 

809.95 

848.62 

91.5 

2953.70 

4519.72 

79.0 

301.50 

303.73 

88.8 

825.98 

866.63 

91.6 

3055.07 

4792.26 

79.2 

306.85 

309.20 

86.9 

842.56 

885.73 

91.7 

3160.32 

5090.87 

79.4 

312.38 

314.82 

87.0 

859.88 

905.41 

91 8 

3209.40 

5418.31 

79.8 

318.03 

320.82 

87.1 

877.38 

925.93 

91.9 

3382.48 

5777.8 4 
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Tabic 1 (contd) 


z 

R 

R' 

92.0 

3499.59 

6174.23 

92.1 

3021.06 

6612.15 

92.2 

3746.20 

7097.56 

92.3 

3875.53 

7638.78 

92.4 

4(8)9.10 

8246.36 

92.5 

4147.07 

8926.70 

92.0 

4289.48 

9692.48 

92.7 

4436.33 

10560.24 

92.8 

4587.50 

11551.26 

92.9 

4742.84 

12684.70 

93.0 

4902.77 

13986.89 


Tabic 2. Maximum refraction model residuals for selected 
P, T, and ranges of Z 


Temperature, 

C 

Maximum refraction model residuals, see 

P - 700 

P - 760 

P - 800 

a. 0 < Z < 85 

-10 

+ 4.59 

+ 4.87 

+ 5.05 

0 

-I 5.25 

+ 5.59 

+ 5.82 

+ 10 

+ 5.83 

+ 6.24 

4 6,51 

+ 20 

4 6.41 

+ 6.88 

+7.18 

4 30 

+ 6.98 

+ 7.50 

+7.83 

b. 85 < Z< 93 

-10 

-15.41 

- 18,36 

-24,30 

0 

13.56 

- 15.03 

— 17.45 

+ 10 

— 1 1.91 

-14.27 

-14.02 

4 20 

-15.16 

-14.77 

— 12.61 

4 30 

- 19.20 

- 16.20 

+ 13.94 


Table 3. Maximum refraction model residuals for selected 
P, T, and ranges of Z: A,, A,, A, - 0 


Temperature, 

Maximum 

refraetion model residuals, see 

"C 

P = 700 

P = 760 

P = 800 

a. 0 < Z < 85 

-10 

+ 6.06 

4 6.38 

+6,57 

0 

+ 5.33 

+ 5.61 

+ 5.78 

+ 10 

+ 4.68 

+ 4.93 

+ 5.08 

+ 20 

+ 4.11 

+ 4.32 

+4.45 

+30 

+ 3,59 

—5.41 

6.06 

b. 85- £ Z < 93' 

-10 

+ 102.05 

-188.07 

278.89 

0 

-130.49 

-251.98 

-338.33 

+ 10 

-196.42 

312.90 

- 395.56 

+20 

-258.65 

370.69 

450.09 

+30 

-317.28 

• 452.40 

501.90 
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Table 4. Surface refractivlty vs integrated rafractivity : A = 0.3224; * = 0.93% 


Case 

MV, 

100 x nd; 
(%) 

( MV,) 

Ax t 100x nd;) 
(%) 

SNW(h)dh 

/ \’D(/« ) dli 
(%) 

A(%) 

4, cm 

1 

4.86 

1.57 

2.27 

+0.70 

+ 1.48 

2 

3.76 

1.21 

2.17 

+0.96 

+2.03 

3 

5.74 

1.85 

1.80 

-0.05 

-0.11 

4 

5.34 

1.72 

1.37 

-0.35 

-0.74 

5 

4.74 

1.53 

1.75 

+0.22 

+0.47 

• 

7.14 

2.30 

2.17 

-0.13 

-0.28 

7 

24.11 

7.77 

8.55 

+0.78 

+ 1.85 

8 

31.72 

10.23 

9.12 

— 1.11 

-2.35 

9 

7.29 

2.35 

4.58 

+2.23 

+4.72 

10 

9.89 

3.19 

2.69 

-0.50 

-1.06 


It 
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Fig. 5. Least squares fit of In R (Garfinkel data) to an 
nth degree polynomial 
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130 



Z ACTUAL, dag 

Fig. 9. Refraction model with and without (pressure correction factor) 
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Appendix A 
Subroutine BEND I 


00101 

1* 

SUBROUTINE REND ( PRESS > TEMp 1 7NITH * P ) 

00101 

2* 

DIMENSION A(2).B<2>,C(2>*E<12>»P<2>*T<2)*Z{2> 

00104 

3* 

PCD = 760.00 

00105 

4* 

P<2) = PRESS 

00106 

5* 

TCI) = 273.00 

00107 

6* 

T ( 2 ) = TEMP 

00110 

7* 

Z C 1 ) 5 91.870 

00111 

A* 

Z<2> = 7NITH 

00112 

9* 

A ( 1 ) = .40816 

00111 

10* 

A ( 2) = 112*30 

00114 

11* 

B ( 1 ) = .12820 

00115 

12* 

B ( 2 ) = 142.88 

00116 

13* 

CC1) = .80000 

00117 

14* 

C (2) = 99.344 

00120 

15* 

E(l) = 46.6?5 

00121 

16* 

E ( 2 ) = 45.375 

00122 

17* 

EC 3) = 4.1572 

00121 

18* 

E ( 4 ) = 1.4468 

00124 

19* 

E(5) = .25391 

00126 

20* 

E ( 6) = 2.2716 

00126 

21* 

E(7) =-1.3465 

00127 

22* 

E (8) =-4.3877 

00130 

23* 

E (9) = 3.1464 

0U131 

24* 

E ( 1 0 ) = 4.5201 

0013? 

25* 

E(ll)=-1.89«2 

00131 

26* 

E ( 12 > = .89000 

00134 

27* 

D3=1.*DELTA(Z»C»Z(2) ) 

00136 

28* 

FP=(P(2)/P(1) )*(1. -DELTA! o»Ar7(?) )/03) 

00136 

29* 

FT=(Td)/T<2> )* d. -DELTA (T*R» 7(2) )/D3) 

00137 

30* 

ll= c z c 2 > -E C 1 ) >/E<2) 

00140 

31* 

X=E(11) 

00141 

32* 

no l i=i»8 

00144 

33* 

1 X=Edl-I )*U*X 

00146 

34* 

R=FT*FP*(EXP(X/D3)-e-(l2) ) 

00147 

35* 

return 

00150 

36* 

END 


00101 

1* 

function delta(a»b»7> 

00103 

2* 

DIMENSION A ( 2 ) »9<2> 

00104 

3* 

DELTA= ( A ( 2 ) - A ( 1 ) )*ExP(Bd)*(7-B(2) ) ) 

00105 

4* 

return 

00106 

5* 

end 
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Appendix B 
Subroutine BEND II 


00101 

1* 

SUBROUT 1 NE BEND < PRESS * TEMP * ZNI TH * R ) 

00103 

2* 

DIMENSION EC 12) 

00104 

3* 

P = 760*00 

00105 

4* 

T = 273.00 

00106 

5* 

E(l) = 46.625 

00107 

6* 

E ( 2) = 45.375 

00110 

7* 

E(3> = 4.1572 

00111 

8* 

E(4) = 1.4468 

0011? 

9* 

E < 5) = .25391 

00113 

10* 

E ( 6) = 2.2716 

00114 

11* 

E ( 7) =-1.3465 

00115 

12* 

E ( 8) =-4.3877 

00116 

13* 

E ( 9) = 3.1484 

00117 

14* 

E ( 10 ) = 4.5201 

00120 

15* 

E ( 11 ) =-1.8982 

00121 

16* 

E( 12) = .89000 

00122 

17* 

FP=PRESS/P 

00123 

IB* 

FT=T/TEMP 

00124 

19* 

U=(ZNITH-E(1))/E<2) 

00125 

20* 

X=E(11) 

00126 

21* 

DO 1 1=1*8 

00131 

22* 

1 X=EU1-I)*U*X 

00133 

23* 

R=FT*FP*<EXP(X)-E(12) ) 

00134 

24* 

return 

00135 

25* 

END 
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Appendix C 
Subroutine SBEND 




00101 

1* 

SUBROUTINE SBEN0CPRESS»TEMP»HUMID»ZNXTH»R) 

00103 

2* 

DIMENSION AC2)*BC2)#CC2)*Efl2)*PC2)*TC2)*ZC2) 

00100 

3* 

PCI) s 760.00 

00109 

0* 

TCI) = 273.00 

00106 

5* 

ZC1) = 91.870 

00107 

6* 

PC2) = PRESS 

00110 . 

7* 

TC2) = TEMP 

00111 

8* 

ZC2) 7 ZNITH 

00112 

9* 

AC1) = .00816 

00113 

10* 

AC2) s 112.30 

00110 

11* 

BC1) = .12820 

00115 

12* 

BC2) = 102.89 

00116 

13* 

CC1) s .80000 

00117 

10* 

CC2) = 99.300 

00120 

15* 

EC1) = 06.625 

00121 

16* 

EC2) = 05.375 

00122 

17* 

EC3) = 0.1572 

00123 

IS* 

ECO) s 1.0068 

00120 

19* 

EC5) = .25391 

00125 

20* 

EC6) s 2.2716 

00126 

21* 

EC7) s-1.3065 

00127 

22* 

ECO) =-0.3877 

00130 

23* 

EC9) = 3.1080 

00131 

20* 

EC10)= 0.5201 

00132 

25* 

EC1D=-1.B982 

00133 

26* 

EC 12) 7 .89000 

00130 

27* 

WO = 7100.0 

00135 

28* 

W1 = 17.109 

00136 

29* 

W2 = 0680.1 

00137 

30* 

W3 = 38.050 

00100 

31* 

D3=1.*0ELTACZ»CZC2) » 

00101 

32* 

FP=CPC2)/PC1))*C1— 0ELTA(P»A»ZC2))/03) 

00102 

33* 

FT=CTCl)/TC2))*Cl.-0ELTACTtBtZC2))/D3) 

00103 

30* 

FW=l*<WO*HUNlD*EXP< CW1*TC2)-W2)/CTC2)-W3) )/CTC2)*PC2) ) > 

00100 

35* 

U=CZC2)-EC1))/EC2) 

00105 

36* 

XsECll) 

00106 

37* 

00 1 1=1*0 

00151 

38* 

1 X=EC11-D*u*x 

00153 

39* 

R=FT*FP*FW* C EXP C X/D3 ) -€ C 12 ) ) 

00150 

00* 

RETURN 

00155 

01* 

END 

00101 

1* 

FUNCTION OEtTACA*B*Z) 

00103 

2* 

DIMENSION AC2)*BC2) 

00100 

3* 

0ELTA=CAC2)-AC1))*EXP(BCD*CZ-8C2))) 

00105 

0* 

RETURN 

00106 

5* 

END 
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Appendix D 
Subroutine XBEND 


00)01 

)• 

SUBROUTINE X*END< PRESS* TCHPtHUHlO*ZN|fH.R) 

0010) 

2* 

dimension e i 1 2 » 

00104 

)• 

P • 740. OU 

001 Oft 

4* 

T • 273.00 

00104 

»• 

C(|) • 44.42ft 

00107 

4* 

C 1 2 ) ■ 4ft.)7S 

00)10 

7 • 

EDI ■ 4 • 1 572 

00))| 

• • 

C C 4 | ■ |.444* 

00)12 

4* 

Oft) • *25)41 

00)1) 

10* 

04) • 2*27)4 

00)11 

))• 

07) ■•1*2445 

0011ft 

|2* 

OB) ■•4*)B77 

00)14 

))• 

04) ■ 3*14*4 

00)17 

14* 

OIOI« 4*5201 

00)20 

16* 

O 1 1 )«-|»*4*2 

00)2| 

14* 

O 1 2 ) ■ **4000 

00)22 

17* 

■0 • 7)00*0 

00)2) 

)•• 

• 1 • 1 7 * ) 44 

00)24 

14* 

•2 ■ 44*4*1 

00) 2ft 

20* 

•) • )• • 450 

00)24 

2 ) • 

FP-PRESS/P 

00)27 

22* 

FW/TCHP 

00))Q 

2)* 

F**l*MO*HUH|OeEXP( («l •TEHP.N2)/ 1 TbHP«*i| )/ ( TEHP*PR£SS) 

00))) 

24* 

U«(2NI TH»0 1 ) 1/02) 

00) )2 

2ft* 

X«dll) 

00))) 

24* 

00 ) 1-1.4 

00) )4 

27* 

1 X»E ( 1 1 • | 1 *U*X 

00)40 

24* 

P.rT*FP*F**(tiP(X |-E( |2) ) 

00)4) 

2f* 

RETURN 

00)42 

)0* 

End 
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